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Abstract 

Based on the tetrahedral Zamolodchikov algebra, we prove the Yang-Baxter equation 
for the R-matrix of 1-D SU (n) Hubbard model. Furthermore, we present generalizations 
of the model. 
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1 Introduction 

The Hubbard model is one of the significant models in the study of strongly correlated 
electronic systems which might reveal an enlightenning role in understanding the mysteries 
of the high-Tc superconductivity. The 1-D Hubbard model also favours a lot of properties of 
integrable models in non-pertubative quantum field theory and mathematical physics. Since 
Lieb and Wu|j^] solved the 1-D Hubbard model by Bethe ansatz method in 1968, based on 
their results (Lieb and Wu's Bethe ansatz equations), many works||2|, ^, ^, |5|, ^, |^, ^ ^, IC, 



11, 12, 13, |l^ have been done extensively to clarify the physical properties of this model. 
Although there were lots of works on the Hubbard model, the integrability was finished until 
1986 by ShastrypI, Olmedilla and Wadati|l|] in both boson and fermion graded versions. 



However, the Yang-Baxter equation for the iZ-matrix of the model obtained by Shastry was 
proved until 1995 by Shiroishi and Wadati in Refs. 32, 33 1 and a generalization of the 



Shastry's bilayer vertex model was also presented in pl|. Moreover, the eigenvalue of the 
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transfer matrix related to the Hubbard model was suggested in Ref. [^] and proved through 
different methods in Refs. |17, 18 1. 

Based on the knowledge of Lie algebra, Maassarani and Mathieu succeeded in constructing 
the Hamiltonian of the SU{n) XX model and showed its integrability|22]. Considered two 
coupled SU (n) XX models, by using Shastry's method, Maassarani constructed the SU{n) 
Hubbard model[^3| and found the related i?-matrix which ensures the integrability of the 1-D 
SU{n) Hubbard model|^^. (It was also proved by Martins for n = 3,4.|25], and by Yue and 
Sasaki|26] for general n in terms of Lax pair formulism.) The exact Solution of the SU{3) 
Hubbard model was given in Ref.[27|. But the Yang-Baxter equation for the given i?-matrix 
was not proved. 

The main purpose of the present paper is to prove the Yang-Baxter equation for the 
i?-matrix of 1-D SU{n) Hubbard model following method suggested in Ref.|31|. In section 
2 we review the model and its integrability. We present the L-operator and the i?-matrix of 
the model and formulate the Yang-Baxter relation. In section 3 we construct the tetrahedral 
Zamolodchikov algebra related to the SU{n) Hubbard model. The Yang-Baxter equation for 
the cooresponding i?-matrix was proved in section 4 and we also present a generalization of 
the model in this section. In section 5 we make some conclusion remarks. 



2 The 1-D SU{n) Hubbard model and its integrability 

The Hamiltonian of the 1-D SU{n) Hubbard model is: 



L n-l 



Un 



2 L 



^ - 2^ Z^{^(T,k^a,k+l + ^a,k^a,k+l + ^T,k^T,k+l + ^T,k^T,k+l) + " Z^^k ^fc ' U 



k=l a=l 



k=l 



where U is the Coulumb coupling constant, and -^"^,(0 = cr,T) is a matrix with an one at 
row a and column (3 and zeros otherwise: 



Im 



The subscripts a and k stand for two different E operators at fc-th site [k = 1, • • • ,L). The 
n X n diagonal matrix C is defined by C = J2a<n ~ We also assume the periodic 
boundary condition: E' 



af3 
k+L 



E' 



In this model, the system has two types of particles named a and r respectively, and each 
particle can occupy (n — 1) possible states. The same type of particles can not appear in 
one site, but two different types of particles can occupy a same site. We denote \n)j as the 
vaccum state of j-th site, \2)j, • • • , |n — l)j as the (n — 1) possible one particle states of 
j-th site. Under the following basis: 

/1\ /0\ /0\ /0\ 







, • • • , |n - l)j 



^0/ 



1 



|n)j 





V 1 / 
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It could be easily proved that: Ef''\n)j = \a) j , E'^°'\n) j = 0,£'j"°|a)j = \n)j,E'^°'\n)j = 0. 
This means that the operators and i?"" can be interpreted as the particle creation and 
annihilation operators respectively. E'^'^ can create a \a)j state particle over the vaccum state 
\n)j of j-th site, and E^"^ annihilate a \a)j state particle to vaccum state of j-th site. 

The SU{n) Hubbard model is constructed by considering two coupled SU{n) XX model, 
so the Hamiltonian consists of two SU{n) XX model with an interation term between 
them. The Hamiltonian of the SU{n) XX model is: 

L n-l 

TT \ ^ \ ^ / TTina Tpan , rjan rpna \ (r,\ 

Hxx - 2^ 2^yEj, i^fc+i + j^fc ^fc+iJ) l^J 

k=la=l 

and the corresponding i?-matrix is: 

R{X) = a(A)[S""(g)S""+ E'^f^^E'^'^] 

a,f3<n 

+h{\) (xS"'^ ® + 

a<n 

+c(A) Y (^''"^"" + ^""^""), (3) 

where x = e^^ and a(A) = cos(A),6(A) = sin(A),c(A) = 1. The functions a(A), 6(A), c(A) 
satisfy the free-fermion relation: a^(A) + 6^ (A) = c^(A). 

The i?-matrix of the SU{n) XX model satisfies regularity property i?(0) = P, unitarity 
condition i?i2(A)i?2i(— A) = cos^(A)/(i and Yang-Baxter equation (YBE): 

-R3l(Al)-R32(A2)-Rl2(A2 — Ai) = i?l2(A2 — Ai )ii32 ( A2)-R3l (Ai ) , (4) 

where P is a permutation operator on the tensor product of two n-dimensional spaces. It is 
easy to verify that it also satisfies a decorated Yang-Baxter equation (DYBE): 

-R3l(Al)i?32(A2)C2i?l2(A2 — Ai) = i?l2(A2 — Ai)C2i?32(A2)i?3l(Ai). (5) 

For the two SU (n) XX model, the i?-matrix without interaction term is given by: 

ii.,(A) = iiJf(A)i?;;)(A), (6) 

here i?-J''(A) and i?jj''(A) denote the ii-matrices of two SU{n) XX models. Since both r\j\x) 
and R'fj\x) satisfy the YBE and DYBE, the product Rij{X) also satsfy the YBE: 

-R3l(Al)-R32(A2)-Rl2(A2 — Ai) = -Rl2(A2 — Ai)i?32(A2)-R3l(Al) (7) 

and DYBE: 

-R3l(Al)-R32(A2)C'2'^''C2^''-Rl2(A2 — Ai) = -Rl2(A2 — Ai )C2'^^ C2^^^32 (A2)-R31 (^1 ) (8) 
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A linear combination of (^) and (^) yields: 

ii3i(Ai)i?32(A2) {aRu{X2 - Ai) + pC^''^ C^^'^ Ri2{X2 + Ai)} 

= {aRi2{X2 - Ai) + (3Ri2{X2 + Ai)C^"^C^")} i?32(A2)i23i(Ai), (9) 

here a and /? are combination coefficients and arbitrary. 

For the SU (n) Hubbard model, the two coupled SU (n) XX model, we look for a solution 
of the Yang-Baxter relation (YBR): 

/:3l(Al)/:32(A2)i??2(Al, A2) = i??2(Al, A2)/:32(A2)/:3l(Al) (10) 

in the form: 

RuiXi, A2) = ai?i2(A2 - Ai) + /?i?i2(A2 + Xi)ct'^ci^\ (11) 
Cij{X) = RijiX) exp{/i(A)cf ^Cf (12) 

Comparing eq.(|9|) with the Yang-Baxter relation (|lO|), we get a relation: 

/i(Ai)/2(A2)i?^2(Ai, A2)/ri(Ai)/2-^(A2) = ai?i2(A2 - Ai) + pd^""^ C^^^ Ri2{X2 + Ai), (13) 
where 

/,(A)=exp{/i(A)c]'^)C7f }. (14) 

Prom (p!3[), we have: 

/3a(A2 + Ai)c(A2 + Ai; 



aa(A2 - Ai)c(A2 - Ai) 
/3b(A2 + Ai)c(A2 + Ai) 
a6(A2 - Ai)c(A2 - Ai) 



tanh(/i(A2) - h{Xi)), 
tanh(/i(A2) + /i(Ai)), (15) 



which give the ratio of a and /? and constraints on h(Xi) and h(X2)- The constraints can be 
written in more explicit form p^ ]: 

sinh(/i(Ai)) ^ sinh(/t(A2)) ^ 
sin(2Ai) sin(2A2) 4 ' ^ ' 



Now we have obtained the i?-matrix of the 1-D SU{n) Hubbard model |24]: 
i?^2(Ai, A2) = 4? (A2 - Ai)4J (A2 - Ai) + ~ M tanh(/i(A2) - /i(Ai)) 

COS(^A2 + Alj 

x4^)(A2 + Ai)4J(A2 + Ai)C^")c^"\ (17) 

which satsify the Yang-Baxter relation (|To|). This i?-matrix depends not only on the diference 
of the spectral parameters A2 — Ai, but also on the sum of the spectral parameters A2 + Ai. 
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This non-additive property allows us to generalize the Hamiltonian of the 1-D SU (n) Hubbard 
model (see section 4). 

The monodromy matrix of the model can be defined as: 

Ta{X) = CLa{X)CL-la{X) ' ' ' Aa(A). (18) 

From the Yang-Baxter relation (|l^) we know the monodromy matrix satisfies the global 
Yang-Baxter relation: 

Ti(Ai)T2(A2)ii5'2(Ai, A2) = i?^2(Ai, A2)T2(A2)Ti(Ai). (19) 

The corresponding transfer matrix is defined by: 

T{\)=tra[Ta{X)] (20) 

and from ( [l9|) it can be easily proved the existence of a commuting family of transfer matrices 

[r(Ai),r(A2)]=0. (21) 

Then the integrability of the model was proved. 

Using the relation /i(0) = and /i'(0) = we can obtain the Hamiltonian of the 1-D 
SU{n) Hubbard model (|): 

H = i-lnT(A)U=o = r-^(0)|^r(A)U=o 

L n-l 

2^[^a,k^a,k+l + ^a,k^a,k+l + ^T,k^r,k+1 + ^T,k^T,k+l) 
k=l a=l 

.2 L 

I 

k=l 

3 Tetrahedral Zamolodchikov algebra 



(22) 



In the above section we have shown the integrability of the 1-D SU{n) Hubbard model. 
It is natrual to expect that the i?-matrix (W% itself satisfy the Yang-Baxter equation(YBE): 



-^3l(A3> Ai)i232(A3, Ai)i?5'2(Ai, Ai) — ii^2(Al> A2)-R32(A3, A2)-R3i(A3, Ai). (23) 



In the SU{2) case, the YBE of the i?-matrix was proved in ref.|31| by using the tetrahedral 



Zemolodchikov algebra (TZA) p4| , 33]. In this section, we construct the TZA related to the 
SU{n) Hubbard model. 

The TZA is defined by the following set of relations: 

r-a ph fC _ \^ nahc ne pd (nA\ 
def 
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where a, 6, •••,/ = 0, 1 and 5f ^ are some scalar coefficients. 
We take >C^^ and >Cj^ as follows: 

^% = Rjki^k - Aj), Cjj^ = Rjki^k + ^j)Ck, (25) 

where Rjk{X) is the i?-matrix of the SU{n) XX model as before. Then we could find the 
following relations which give the TZA (p^): 

nO rO rO rO rO rO rO rl rl rl rl rO /'9R'\ 

''-12''-32''-31 — ''-31''-32''-12 5 ''-12''-32''-31 — ''-31''-32''-12' l^^j 

rl rO rl rl rl rO rl rl rO rl (07\ 

'-'12'-'m'-?,l — ''-31''-32''-12i ''-12''-32''-31 " ''-31''-32''-12 ' l^V 

rl rl rl _ nlU rl rO rO , qIU rl rO rl , qIU rO rO rl f^Q\ 

''-12''-32''-31 — '-'001''-31''-32''-12 "r '-'010''-31''-32''-12 + '-'lOO''-31''-32''-125 K^") 

rO rO rl _ qOOI /-l rl rl . qOOI rO rO rl , qOOI rO rl rO /r,o\ 

''-12''-32''-31 — '-'lll''-31''-32''-12 '-'l00''-31''-32''-12 + '-'010''-31''-32''-125 l^^j 

rO rl rO _ qOIOz-I rl rl , qOIO rO rO rl , qOIO rl rO rO fon\ 

''-12''-32''-31 — '-'lll''-31''-32''-12 "r '-'l00''-31''-32''-12 '-'OOl ''-31''-32''-12 1 K'^^J 

rl rO rO _ c-lOOz-l rl rl , clOO /'O rl rO , qIOO rl rO rO /o-\\ 

''-12''-32''-31 — '-'lll''-31''-32''-12 "r '-'OlO ''-31''-32''-12 "r '-'OOl ''-31''-32''-12 ' K'^^J 

where the coefficients S^^j are given by 

sinO 

o^ 

S 

S' 

s: 



.111 


sin(A2 + Ai) cos(A2 + A3) 


olll 

'-'010 — 


sin(A2 + Ai) sin(Ai + A3) 


001 — 


cos(A2 — Ai) sin(A2 - A3) ' 


cos(A2 — Ai) cos(Ai — A3) 


.111 


sin(Ai + A3) cos(A2 + A3) 


cOOl 

'-'ill - 


sin(A2 - Ai) cos(A2 - A3) 


100 — 


cos(Ai - A3) sin(A2 - A3) ' 


cos(A2 + Ai) sin(A2 + A3) ' 


.001 


sin(A2 - Ai) sin(Ai + A3) 


cOOl 

'-'010 — 


sin(Ai + A3) cos(A2 - A3) 


100 — 


cos(A2 + Ai) cos(Ai - A3) ' 


cos(Ai - A3) cos(A2 + A3) ' 


.010 


sin(A2 — Ai) sin(Ai - A3) 


cOlO 

'-'100 — 


sin(A2 - Ai) cos(A2 + A3) 


111 — 


cos(A2 + Ai) cos(Ai + A3) ' 


cos(A2 + Ai) sin(A2 - A3) ' 


.010 


sin(Ai - A3) cos(A2 + A3) 


qIOO 

'-'ill — 


sin(Ai - A3) cos(A2 - A3) 


001 — 


cos(Ai + A3) sin(A2 - A3) ' 


cos(Ai + A3) sin(A2 + A3) ' 


.100 


sin(A2 + Ai) cos(A2 - A3) 


c-100 
'-'OOl — 


sin(A2 + Ai) sin(Ai - A3) 


010 — 


cos(A2 - Ai) sin(A2 + A3) ' 


cos(A2 - Ai) cos(Ai + A3) 



r.100 _ -r y^L) ^^oyyvz - ^V3; ^-loo _ oiiiv/v2 -r y^ij Diiiy/vi — /V3; , , 

" cos(A2 - Ai) sin(A2 + A3) ' " cos(A2 - Ai) cos(Ai + A3) ' ^ ' 

Eg. ([2^) and eq.(p7|) are equivalent to the YBE (0) and DYBE (^) respectively. In this sense, 
the TZA ( [2^ ) can be regarded as a generalization of the YBE and DYBE. 

It is important to notice that the products £f2'^32^3i ^'^^ linearly independent as 
operators acting on Vi V2 V3 and they satisfy the following relations: 

'^12-^32'^31 = ^0'Cl2'^32'^31 + yo-^12'^32'^31 + ^0'Cl2'^32'^31 ) (33) 
'^12'^32'^31 = ^l'^12'^32'^31 + yi'^12'^32'^31 + ^l'^12-^32'^31) (34) 



with 



_ cos(Ai - A3)sin(A2 - A3) _ cos(A2 - Ai) cos(Ai - A3) 

cos(Ai + A3) sin(A2 + A3) ' cos(A2 + Ai) cos(Ai + A3) ' 
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_ cos(A2 - Ai)sin(A2 - A3) _ cos(Ai + A3) sm(A2 + A3) 
cos(A2 + Ai) sm(A2 + A3) ' cos(Ai - A3) sm(A2 - A3) ' 

_ cos(A2 + Ai) cos(Ai + A3) _ sin(A2 + A3) cos(A2 + Ai) 
cos(A2 — Ai) cos(Ai — A3) ' sin(A2 — A3) cos(A2 — Ai) ' 

From these relations we know that the hnear space spanned by the products £^2^32^31 
is 6-dimensional. 



4 The Yang-Baxter equation for the i?-matrix of the SU{n) 
Hubbard model 

In this section we prove the Yang-Baxter equation for the i2-matrix of the 1-D SU (n) 
Hubbard model (H). 

Taking into account the form of the i?-matrix (0), we look for a solution of the YBE 
(|2^ in the following form: 

R%{X„ Afc) = ii]?(Afc - A,)i?J^ (A, - A,) + a,kR{<r)^ki^k + A,)cf ^ijj^ (A, + X,)ci^'^ 

= <r^<r + a,.4r^4r, (38) 



where C^i^ and have been defined in (|^). If ajk = 0, the i?-matrix satisfies the YBE (|2^ 
in a trivial way. Now we look for a non-trivial solution. Subsitituting the expression 
into the Yang-Baxter equation (|2^) , by means of the tetralhedral Zamolodchikov algebra and 
relations ( [3^ ) and (^4|), we could find that ajk must satisfy the following conditions: 



If we take 



ai2 sin 2(Ai -|- A2) + a^i sin 2(Ai -|- A3) = 032 sin 2(A2 -|- A3) 

= — sin2(A3 - Ai) — sin2(A2 - Ai). (39) 
"31 "12 

= 7 tanh(/i(Afc) - h{Xj)), (40) 

cos(Afe + Xj) 



and impose the constrains 

sinh(/i(Aj) _ n^U 



(i = 1,2,3), (41) 



sin(2Aj) 

then the conditions ( |39|) are satisfied. This proves the Yang-Baxter equaiton for the i?-matrix 
of the 1-D SU (n) Hubbard model. 

Besides the YBE (^), the i?-matrix ( p^ has the following properties: 

1 

coshJhiX))' 

R%iXo,Xo) = V,k, (43) 
R^iXj, Xk)R%{Xk, Xj) = p{X,,Xk)Id, (44) 
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where 

p{\j, Afc) = cos^(Afc - Xj) |cos^(Aa; - \j) - tanh^(/i(AA:) - /i(Aj))| , (45) 
and the permutation operator is defined as 

T'^k = vf^vfl (46) 
The Yang-Baxter equation ( p3| ) imphes a more general inhomogeneous model as: 

Ta{\ {A,}) = i?L(A, \n)RI-U\ Al_i) • • • <(A, Ai), (47) 
where Aj, (j = 1, 2, • • • , L) are the inhomogeneous parameters obeying the constraints 

0^1.2,....L), ,48) 

sm(zAjj 4 

From the Yang-Baxter equation (p3|), we can obtain the global Yang-Baxter relation: 

Ti(A, {A,})T2(^, {A,})i?^2(A, /u) = RU\ {A,})Ti(A, {A,}), (49) 

which leads to the commutativity 

[T(A,{A,}),r(/x,{A,})]=0, (50) 

where r(A, {Aj}) is the transfer matrix of the model 

T(A,{A,})=tr„T„(A,{A,}). (51) 

The corresponding Hamiltonian is defined as the logarithmic derivative of the transfter matrix 
under all inhomogeneous parameters Xj = Aq, (j = 1, 2, • • • , L): 

^Ao = ^^'^'^(^'{^J = ^o})|a=Ao = T-~^{Xo,{Xj = Xo})-^T{X,{Xj = Ao}) 
L 

j=l a<n 

+ 4cosh(2/i(Ao)) ^^^^ 

where 

= cos{2Xo){-EY''E^^,+ J2 E;^"E^!l,) + sm{2Xo)J2{E;^Ef^, + E^^E:i^,) 

a,l3<n OL<n 

+ E (-^^r^iTi + ^r^iTi ) • (53) 

The arbitrariness of the parameter Aq comes from the non-additive property of the spectral 
parameters. If we take Aq = 0, this new Hamiltonian reduces to (|^). 

Thus, we have obtained a new 1-D SU{n) Hubbard by the Yang-Baxter equation of the 
i?-matrix (23). 
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5 Conclusions 



In this paper we have proved the i?-matrix of the 1-D SU{n) Hubbard model satisfying 
the Yang-Baxter equation. We notice that the tetrahedral Zamolodchikov algebra play an 
essential role in the proof. 

In most lattice systems, the existence of i?-matrix ensures the integr ability and the R- 
matrix is isomorphic to the L-operator. Thus, the Yang-Baxter equation is a consequence 
of the Yang-Baxter relation R12L1L2 = -L2^i-Ri2- But, for Hubbard model, the situation 
is quite different. The i?-matrix can not be obtained from L-operator, even if we limit to 
SU{2) Hubbard model[15, |^]. The ii-matrix satisfying the Yang-Baxter equation together 
with L-operator constitute the complete proof of the integrability. 

For SU (n) Hubbard model, the i?-matrix does not isomorphic to the L-operator. This 
provides a method to construct a new kind of integrable system by considering the i?-matrix 
as a L-operator (fundamental representation of same algebra). The general representation 
can be obtained by fusing the multi fundamental rep. (i?-matrix). Therefore, one can get 
the full representation of the algebra in principle. 

In the present paper we have derived out a new Hamitonian ( |52|) from the ii!-matrix. The 
last term introduces a new kind of interaction. This Hamiltonian is quite different from the 
original one (||). It also renders a question how to find the eigenvalue of this Hamiltonian 
(jS^. We will consider it late. 

In the derivation of eq. (|52|), we have assumed all the parameters \j are same Aq. This 
is not necessary. The different choice will give out different Hamiltinian. On the other hand, 
one can consider the i?-matrix as a Boltzmann weight in statistical mechanics. This provides 
a inhomegeneous lattice statistical medel. The partition function could be derived out in a 
similar way. 
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